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' Abstract The clothing procedure, put forward in quan- 
] turn field theory (QFT) by Greenberg and Schweber, is 

■ applied for the description of nucleon-nucleon (TV — N) 
scattering. We consider pseudoscalar (tt and rj), vec- 
tor (p and uj) and scalar {S and a) meson fields in- 

, teracting with 1/2 spin (N and N) fermion ones via 
' the Yukawa-type couplings to introduce trial interac- 
. tions between "bare" particles. The subsequent uni- 
] tary clothing transformations (UCTs) are found to ex- 

■ press the total Hamiltonian through new interaction 
operators that refer to particles with physical (observ- 

■ able) properties, the so-called clothed particles. In this 
, work, we are focused upon the Hermitian and energy- 
' independent operators for the clothed nucleons, being 
, built up in the second order in the coupling constants. 
' The corresponding analytic expressions in momentum 

■ space are compared with the separate meson contribu- 
' tions to the one-boson-exchange potentials in the me- 
son theory of nuclear forces. In order to evaluate the T 

, matrix of the N — N scattering we have used an equiva- 
' lence theorem that enables us to operate in the clothed 
particle representation (CPR) instead of the bare par- 
ticle representation (BPR) with its large amount of vir- 
tual processes. We have derived the Lippmann-Schwin- 
' ger type equation for the CPR elements of the T— ma- 
trix for a given collision energy in the two-nucleon sec- 
tor of the Hilbert space H of hadronic states. 
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1 Introductory remarks and some recollections 

We know that there are a number of high precision, 
boson-exchange models of the two nucleon force Vnn, 
such as Paris [T], Bonn [5], Nijmegen [3], Argonne [1], 
CD Bonn [5] potentials and a modern family of covari- 
ant one-boson-exchange (OBE) ones [S].Note also suc- 
cessful treatments based on chiral effective field theory 
[2], [8], for a review see [9]. 

In this paper, we would like to draw attention to 
the first application of unitary clothing transformations 
(UCTs) [TnUTT] in describing the nucleon-nucleon scat- 
tering. Recall that such transformations being aimed 
at the inclusion of the so-called cloud or persistent ef- 
fects, make it possible the transition from the bare- 
particle representation (BPR) to the clothed-particle 
representation (CPR) in the Hilbert space Ti of meson- 
nucleon states. In this way, a large amount of virtual 
processes induced with the meson absorption/emission, 
the A^A^— pair annihilation/production and other cloud 
effects can be accumulated in the creation (destruction) 
operators ac for the "clothed" (physical) mesons and 
nucleons. Such a bootstrap reflects the most significant 
distinction between the concepts of clothed and bare 
particles. 

In the course of the clothing procedure all the gen- 
erators of the Poincare group get one and the same 
sparse structure on V. [TU] . Here we will focus upon one 
of them, viz., the total Hamiltonian 



H = Hpia) + Hi{a) = H{a) 
with 



(1) 
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Hi{a) ~ V(a) + mass and vertex counterterms, (2) 

where free part Hp{a) ^ a^a belongs to the class [1.1], 
if one uses the terminology adopted in [10], and inter- 
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action V{a) is a function of creation (destruction) op- 
erators a^(a) in the BPR, i.e., referred to bare particles 
with physical masses where they have been intro- 
duced via the mass-changing Bogoliubov-type UTs. To 
be more definite, let us consider fermions (nucleons and 
antinucleons) and bosons (tt— , rj—, p—, oj-mesons, etc.) 
interacting via the Yukawa- type couplings for scalar (s) , 
pseudoscalar (ps) and vector (v) mesons. Then, as seen 
from Appendix A, V{a) = Vs -\- Vps + K with 



(3) 



Vs ^ gs / (ixV'(x)'0(x)(p5(x) 



Vps ^igps / (ixV'(x)75V'(x)(^ps(x) 



(4) 



projection onto the quantization axis ( the particle mo- 
mentum ) for the so-called canonical ( helicity-state ) 
basis ( see, e.g.. Chapter 4 of Ref. [M] ). 

In the context, we have tried to draw parallels with 
that field-theoretic background which has been em- 
ployed in boson-exchange models. First of all, we imply 
the approach by the Bonn group _2, 5 , where, following 
the idea by Schiitte [T5], the authors started from the 
total Hamiltonian (in our notations). 



H = Hpia) + V{a) 

with the boson-nucleon interaction 

V{a) r^ab^ + H.c. 



(8) 



(9) 



K = / dx <{ 5v?/'(x)7mV'(x)¥'(:(x) 

/v 



+ ^t/-(x)a^.V(x)^r(x) 



dx <{ — i^V'(x)7o'/'(x)V'(x)7oV'(x) 



+ ^V5(x)ao.V(x)^(x)ao.V(x) \ , (5) 



where (/^(^''(x) = d^(p'^{x) — 9''(^(f(x) the tensor of the 
vector field included. The mass (vertex) counterterms 
are given by Eqs. (32)-(33) of Ref. [TT] (the difference 
Vo{a) - V{a) where a primary interaction Vo{a) is de- 
rived from V{a) replacing the "physical" coupling con- 
stants by "bare" ones). 

The corresponding set a involves operators a^(a) 
for the bosons, 6^(6) for the nucleons and d^{d) for the 
antinucleons. Following a common practice, they ap- 
pear in the standard Fourier expansions of the boson 
fields ipb and the fcrmion field ip, though for our pur- 
poses the use of such a representation is not compulsory 
(see, e.g.. Chapter 3 of the monograph [H]). In any case 
we have the free pion and fcrmion parts 

Hpia) = J dka;ka^(k)a(k) 

dp Yl (P' ^)^(P' + '^^ (P' ^)^(P' ^)] (6) 



For brevity, the contributions from other mesons are 
omitted. 

Unlike Eq.Q the antinucleonic degrees of freedom 
were disregarded in Ref. [2]. There (see also [T4|) the 
transition matrix T{z) (in the two-nuclcon space) was 
considered in the framework of the three-dimensional 
perturbation theory when handling the integral equa- 
tion, 

TNNiz) = Vnn{z) + VNN{z)iz - HNr^TNNiz), (10) 

where the energy-dependent " quasipotential" Vnn{z) 
approximates a sum of all relevant non-iterative dia- 
grams, Hjf the nucleon contribution to Hp. Such a po- 
tential has "the unpleasant feature of being energy- 
dependent. This complicates applications to nuclear 
structure physics considerably" (quoted from p. 40 in 
[2]). Therefore, further simplifications are welcome (see, 
e.g., Refs. m, [S] ). 

Along with our derivation of a Lippmann-Schwinger 
(LS) equation for the T matrix of the N — N scattering, 
we will demonstrate its solutions to be compared with 
those by the Bonn group. 



2 Analytic expressions for the quasipotentials 
in momentum space 

As shown in [10], after eliminating the so-called bad 
term^ from V{a) the primary Hamiltonian H(a) can 
be represented in the form. 



and the primary trilinear interaction 
V{a) - ab^b + abW' + adb + add^ + H.c. (7) 

with the three-legs vertices. Here cjk = ^/rnl + ^(") = ^pi^^c) + Kj{a^) = K{ac) 

{Ep = \J vn? + p^) represents the pion (nucleon) en- 
ergy with physical mass mii{m) while /j, denotes the 
fermion polarization index. It may be the particle spin 



(11) 



^ By definition, they prevent the bare vacuum 17o {i|^o) 
b\Qo) = ••■ = 0) and the bare one— particle states \lbare) 
at|i7g) Qi'iIQq), . . .) to be -ff eigenstates. 
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The free part of the new decomposition is determined 

by 

KFiac) = J dkuji^al{k)ac{k) 

dp Ep [blip, f^)bc{p, Ai) + 4(P. m)4(p, /i)] 



(12) 

while Kj contains only interactions responsible for phys- 
ical processes, these quasipotentials between the clothed 
particles, e.g., 

K\^\ac) = K{NN NN) + K{NN NN) 
+ K{NN NN) + K{hN hN) + K{hN hN) 

+ K{hh' ^ NN) + K{NN bb') (13) 

A key point of the clothing procedure developed in 
|10| is to fulfill the following requirements: 

i) The physical vacuum (the H lowest eigenstate) 
must coincide with a new no-particle state i7, i.e., the 
state that obeys the equations 

ac{k) \Q) = bc{p,fi) \Q) = dc{p,^l) \n) =0, V k, p, ^ 

(14) 

{f2\n) ^ 1. 

ii) New one-clothcd-particlc states |k)c = aj(k)i7 etc. 
are the eigenvectors both oi Kp and K, 

K{ac)\k), = KF{a,)\k),^LOk\k), (15) 

Kiiac)\k), = (16) 

iii) The spectrum of indices that enumerate the new 
operators must be the same as that for the bare ones . 

iv) The new operators ac satisfy the same commu- 
tation rules as do their bare counterparts a, since the 
both sets are connected to each other via the similarity 
transformation 



(17) 



with a unitary operator W to be obtained as in pCU] . 

It is important to realize that operator K{ac) is the 
same Hamiltonian H{a). Accordingly [10,11] the N- 
N interaction operator in the CPR has the following 
structure: 

K{NN NN) = KbiNN NN), 



Kb{NN ^ NN) = [Y '^Pi ^P2 

xVbil',2';l,2)bUl')bU2')bcil)bc{2), (18) 



where the symbol ^ denotes the summation over the 

nucleon spin projections, 1 — {pi,/ii}, etc. 

For our evaluations of the c-number matrices Vf, we 
have employed some experience from Refs. |10|[TT| to 
get in the second order in the coupling constants (see 
Appendix A) 



^,(1',2';1,2) = 



1 



(27r)3 ^Ep.Ep,Ep,Ep, 



'5(p'i+pi-Pi-P2)t'6(l',2';l,2), (19) 



^^s(l',2';l,2) 
= -4u{p[)uip^) 



{pi — p[)'^ — 



u{p'^)uip2), (20) 



vpsil',2';l,2) 
= %^"(pi)75M(Pi) 



/^2_^2 w(p'2)75^^(P2), 



(21) 



t;v(l',2';l,2)- 



2 [p'l -Pi)^ - m'^ 



/v / / . 

2^^P'^ 


Pi).'! 


u{pi) 


{p' 

2m 


+ P2)'' 


|"(P2) 


fv 1 1 


+ Pl)u 


|w(Pl) 



X^2(p^)^{(p'i+j5(,-pi-p2)7'' 



(p'l +P2 ~Pl -P2)''}u{p2) 



(22) 



where mi, the mass of the clothed boson (its physical 
value) and q — (7^7^. In the framework of the isospin 
formalism one needs to add the factor r(l)T(2) in the 
corresponding expressions. 

One should stress that in the course of our deriva- 
tions the Feynman propagator 

[{p, ~ p[)' - mir' 

arises from adding the noncovariant propagators 
and 

[2wk {Ep,^ - Ep^ - u;i,)f\ 
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Fig. 1 The typical OFPT diagrams with the intermediate boson 
(dashed lines) on its mass shell. 



Such a feature of the UCTs method allows us to use the 
graphic language of the old-fashioned perturbation the- 
ory (OFPT) (see, e.g., Chapter 13 in Schweber's book 
[17]) when addressing the graphs in Fig.l. As noted in 
|llj the graphs within our approach should not be in- 
terpreted as the two time-ordered Feynmam diagrams. 
Indeed, all events in the S picture used here are related 
to the same instant t = 0. Being aware of this, the line 
directions in Fig.l are given with the sole scope to dis- 
criminate between nucleons and antinucleons. The lat- 
ter will inevitably appear in higher orders in coupling 
constants and for other physical processes (e.g., the ttN 
scattering) as it has been demonstrated in Ref. [TT] . 

Further, for each boson included the corresponding 
relativistic and properly symmetrized N — N interac- 
tion, the kernel of integral equations for the N — N 
bound and scattering states, is determined by 

{bl{p[)bM)n\ K,{NN ^ NN) |6t(pi)6t(p2)l2) 

= l/*'-(l',2';l,2)-Fr'(l',2';l,2) (23) 

where we have separated the so-called direct 

Vf'\l', 2'; 1, 2) = -H(l', 2'; 1, 2) - H(2', l'; 2, 1) (24) 

and exchange 

(25) 



yr^(l',2';l,2) = H'-(2',l';l,2) 



terms. For example, the one-pion-exchange contribu- 
tion can be divided into the two parts: 



V;*''(l',2';l,2) 
X 5 (p'l + P2 - 



9l 



1 



(2^)3 ^/E^TE^TE^, 

Pi - P2)w(pi)75w(Pl)"(P2)75M(P2) 
1 



1 



2 I(pi-Pi) 

and 

Fr^(l',2';l,2)-- 
X 5 (p'l + P2 - Pi - 

1 r 1 



iP2 -P2? - ml 



(26) 



9l 



(2^)3 ^/E^TE^TE^, 

- P2) w(Pi)75u(p2)u(P2)75w(Pl) 



{P2 - p'lY - ml {pi - -p'^f - ml 



(27) 



to be depicted in Fig. 2, where the dashed lines corre- 
spond to the following Feynman-like "propagators": 



1 



1 



on the left panel and 



ml {V2-V'2f 



1 



{Vi - P'l)^ " ml {pi - p'^f - ml 



on the right panel. Other distinctive features of the re- 
sult ((23)) have been discussed in [TOl[TT| . Note also that 
expressions p5)) - (l?7l) determine the one-pion-exchange 
part of one-boson-exchange interaction derived via the 
Okubo transformation method in [18 ( cf. [T^ITTj ) tak- 
ing into account the pion and heavier-meson exchanges. 



Direct p^_^ Exchange p j 

Fig. 2 The Feynman-like diagrams for the direct and exchange 
contributions in the r.h.s. of Eq.(23). 



3 The field theoretic description of the elastic 
N—N scattering 

3.1 The T-matrix in the CPR 

Usually in nonrclativistic quantum mechanics (NQM) 
the LS equation for the T operator 

T{E + iQ) = VNN + VNN{E + i{)~h^)-^T{E + iQ) (28) 

with a given kernel Vnn is a starting point in evalu- 
ating the N-N phase shifts. All operators in Eq. (pS)) 
(including the sum /ip of the nucleon kinetic energies) 
act onto the subspace of two-nucleon states, and re- 
main confined in this subspace (the particle number is 
conserved within nonrclativistic approach) . The matrix 
elements {N'N'\T{E + iO)\NN), form the correspond- 
ing T— matrix and, at the collision energy 



E — El + E2 — El 



E. 



27 



the on-energy-shell elements can be expressed through 
the phase shifts and mixing parameters (see below). 
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In relativistic QFT the situation is completely dif- 
ferent. Though one can formally introduce a field oper- 
ator T that meets the equation 

T{E + iO) = Hi + Hi{E + Hf)~^T{E + iO) (29) 

the field interaction as a rule, does not conserve 
the particle number, being the spring of particle cre- 
ation and destruction. The feature makes the problem 
of finding the N-N scattering matrix much more com- 
plicated than in the framework of nonrelativistic ap- 
proach, since now the T matrix enters an infinite set of 
coupled integral equations. 

Such a general field-theoretic consideration can be 
simplified with the help of an equivalence theorem [211 
I20j according to which the S matrix elements in the 
Dirac (D) picture, viz.. 



Sf, = {aK..f2^\S{a)\a\..Qa) 

are equal to the corresponding elements 



S'j,^{al..Q\S{a,)\al...Q) 



(30) 



(31) 



of the S matrix in the CPR. We say "corresponding" 
keeping in mind the requirement Hi) of Sec. 2. The S 
operators in Eqs. (|30p - (PT|) are determined by the time 
evolution from the distant past to the distant future, 
respectively, for the two decompositions 



and 



H = H{a) = Hf + Hi 



H = K{ac) ^ Kf + Ki 



Note that the equality Sji = Sj^ in question becomes 
possible owing to certain isomorphism between the ac 
algebra and the a algebra once the UCTs Woit) = 
exp{iKFt)W exp{~iKFt) obey the condition 

W^d(±oo) = 1 (32) 

The T operator in the CPR satisfies the equation 

T,iotk{E + lO) = Kj 

+ Ki{E + iQ-KF)-^T,ioth{E + iQ) (33) 

and the matrix 

T/, = {f-b\T{E + iQ)\i-b) 

= {f;c\T,ioth{E + iO)\i;c) = T^,, (34) 

where |; 6) ( |; c) ) are the Hf ( Kf ) eigenvectors, may 
be evaluated relying upon properties of the new inter- 
action Ki(ac). The latter has nonzero matrix elements 
only between the clothed-particle (physical) states. Such 



a restriction helps us to facilitate the further consider- 
ation compared to the BPR with its large amount of 
virtual transitions (cf., our discussion in Sect.l ). 

If in Eq.® we approximate Ki by i^f ^ (see Ea.dTSt'). 
then initial task of evaluating the BPR matrix elements 
{N'N'\T{E + iO)\NN} can be reduced to solving the 
equation 

(l',2'|TjvAr(^)|l,2) = {1',2'\Knn\1,2) 

+ (1', 2'\KNNiE + iO- Kf)-^Tnn{E)\1, 2) (35) 

with Knn = K{NN NN). 

Actually, let us employ the relation, 

GF{z)T,ioth{z) = G{z)Ki (36) 

with the two resolvents 

GF{z)^{z-KFy^ G{z) = {z-K)-^ 

In its turn, Ea. p6p can be rewritten as 

oo 

iGF{z)Taoth{z)^ J dT^^'-''^^Ki (37) 



(38) 



To the approximation in question, where 
X ~ if J. + =K2 + Ki 
with 

K2 = K2{ferm) + K2{mes), 

K2{ferm) b\bc + d\dc, K2{mes) ~ ajoc 

and 



K4 - blb%b, + dldld.d, + bldlbdc 

+ a\b],acbc + al.d\.acdc + d],d\acac + i?.c., 

we have 

KiblbWQ) ^ Kiblbl\Q) = K^blbl\Q), 
= blblbcbc + H.c. 



so 



KiblbWn) ~ e''^e-'^^'+^'^^K^blbl\Q). 



But 



l-i{K2 + Ki)T 



--{K2+K,fT^ 



K^blbl\Q) 



with 



{K2 + K,)K^blbl\Q) = {K^ + K^)K^blbl\n), 
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and so on. Here blbl. 

Thus, to the approximation (j38p 

GF{zmz)blbl\Q) = (z - ~ K^y'K^blbllQ) 

or 



TNN{z)\h2) = K^\1,2)+K^{z~K^)-'Tnn{z)\1,2). 

(39) 

Sometimes it is convenient to use the notation |1,2) = 
blbl\f2) for any two nucleon state. Equation ([55)) fol- 
lows from ([Ml) if we take into account the completeness 
condition 

^\NN}{NN\ = 2 (40) 

NN 

and put z — _E + ioH Here the symbol J] means the 

NN 

summation over nucleon polarizations and the integra- 
tion over nucleon momenta. 



3.2 The i?-matrix equation and its 
angular-momentum decomposition. The phase-shift 
relations 

For practical applications, in order to get rid of some 
discomfort in handling the singularity of the resolvent 
{E+iO — K2)~'^, one prefers to work with the i?~matrix 
which is related to the T-matrix by the Heitler equation 
( see, e.g., Sect. 6 of Chapter V in the monograph [H]): 

T{E) = R{E) - i'KR{E)6{E - K^)T{E) (41) 

Thus, for R{E) in our case we obtain 

Rnn{E)\1,2) =Knn\1.2) 



K 



NN 



E-K. 



N 



Rnn{E)\1,2), (42) 



where P denotes the principal value {p.v.) to be applied 
when the integration over the continuous spectrum of 
the operator 

^2=J2j "^'P^P^c (P, f^) be (P, fi) 



is performed. 

^ Henceforth, the infinitesimal shift +iO will be omitted in 
T{E + iO) 



After this, one can write 
(l'2'|i?(£;)|12) = (1'2'|A-Ar^|12) 



t 

34 



;,,.|A-„„|34)<|Mae ,43, 



with R{E) = R{E)/2 and Knn = Knn/2, where the 
operation^ involves the p.v. integration. 

34 

Certainly, the integral equation (j43p has much in 
common with the two-body i?-matrix equation in NQM. 
It is true, unlike the latter, in our case the center-of- 
mass motion is not separated from the internal motion 
that is typical of relativistic theories of interacting par- 
ticles. The kernel of Eq.(|131) is 

(l'2'l Knn |12) = 5 (p'l + p^, - pi - P2) (l'2'l V |12) 
= 6 (p'l + P2 - Pi - P2) 

X (p'lM'lT{,P2M2'r2 l^^l PlMin,P2M2T-2) 

that provides the total momentum conservation in ev- 
ery intermediate state. The subsequent calculations are 
essentially simplified in the center-of~mass system (c.m.s), 
in which we will employ the notations 

|PMlM2,TiT2) = |PMlTl,-PM2T-2) , 



|p ^J■l^J■2,nT2) 

and 



IpVi^I, -pV2^2) 



\^^J'3^J■i,T3T^) = \qn3T3, -qM4T4) , 

respectively, for the initial, final and intermediate states. 
Here the quantum numbers h{t) are the individual spin 
(isospin) projections. 

Using these notations Ea.(|i5|) in the c.m.s. can be 
written as 

{p'fJ,'l^i2,T[T!i \R{E)\pfilfl2,TlT2) 

= (p'MlM2''ri'^2 P/^lM2,Tir2) 

+ J dq{p'n[fi'2,T[T^ \ V\qfi3H4,T3T4) 

(q^3^4,T3T4 |-R(-E)| p^i^2,Tir2) 



E-2Ec 



(44) 



Accordingly Eq. 



(i'2'|y |i2) = - 



2(2^)3 £;p,Sp 
J2[vti^', 2'; 1, 2) - vra', 2'; 2, 1)] (45) 



with 



z;r(l', 2'; 1, 2) = vt{l' , 2'; 1, 2) -f- «,(2', 1'; 2, 1) (46) 
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and 

«r(l',2';2,l) = «*^(2',l';l,2), 

where the separate boson contributions are determined 
by Eqs. (HOll-dlSI) with Pi = p = -p2 and p'^ = p' = 

Following a common practice we are interested in 
the angular -momentum decomposition of Eq.(|44]) as- 
suming a nonrelativistic analog of relativistic partial- 
wave expansions (see [33] and refs. therein) for two- 
particle states. For example, the clothed two-nucleon 
state (the so-called two-nucleon plane wave) can be 
represented as 

|PAiiM2,TiT2) = (^^fii^H2 \SMs^ \TMt^ 

{ImiSMs I JMj ) (p/p) \pJ{lS)M.j, TMt) , (47) 

where, as shown in Appendix B, J, S and T are, respec- 
tively, total angular momentum, spin and isospin of the 
NN pair. All necessary summations over dummy quan- 
tum numbers are implied. 

Of course, when deriving expansions like (j47|) (see 
Appendix B) one needs to take into account a distinc- 
tive feature of any QFT associated with particle pro- 
duction and annihilation. The use of such expansions 
gives rise to the well known JST representation (see, 
e.g., [22]), in which 

{p'J'{l'S')M'j, T'M'j, \V\ pJ{lS)Mj, TMt) 

= y/f^ip' ,P)h'j6M'jMjSs's6T'T6M^MT^ (48) 

{p'J'{l'S')M'j, T'Mlr \R{E) \ pJ{lS)Mj, TMt) 

= Rf,f'^{p',p; E)5j,j5m',m.,5s's5t't5m^Mt , (49) 

so Ea. (|44)) reduces to the set of simple integral equa- 
tions, 



if we are interested in the i?— matrix on the energy shell 
one has inevitably to go out beyond it. 

In addition, one should note that in view of the 
charge independence assumed in this work one has to 
solve two separate equations for isospin values T — 
and r = 1. Regarding the T— dependence of the corre- 
sponding solutions in details one needs to keep in mind 
the typical factor [1 — (—1)'+'^+-^] that is discussed in 
Appendix C (below Eq.(C.24)). It means such a selec- 
tion of the partial R— equations, where the number 
1 + S + T must be odd. 

Properties and are resulted from rotational 
invariance, parity and isospin conservation in combina- 
tion with the antisymmetry requirement for two-nucleon 
states. At the point, let us recall that the clothing pro- 
cedure violates no one of these symmetries (details can 
be found in survey |10j). In particular, it implies that 

('2) 

[3F{ac),Ki(ac)\= 3 F{ac), K] ' {uc) 

^[iferm,K{NN NN)]^Q, (52) 



[y(ae),X/(ac)] 



and 



y{a,),Kf\a,) 

[Vferm, K (N N ^ NN)] = (53) 



[JF(ae),^(ae)] = [a'(ac),^(ae)] =0, 



(54) 



where Jpioc) — J the operator of total angular mo- 
mentum, CP(ac) = J' the operator of space inversion. 
As in Appendix B, we distinguish their fermionic parts 
J/erm and Tferm determined by Eqs. and (|114p . 

Further, apart from the transition with I = V = J 
and l{l') = J ± 1 all different ones are forbidden owing 
to the parity conservation. Therefore, for each JST— 
channel there are only six independent elements H: 



g2 dq 



R 



JS=OT 
J J 



ip',p), R] 



fr'^iP^p), Rj^Y7.Ap',P), 



2{Ep — Eq) 



Wi^^{p',q)Rf4^{q,p) (50) Rr-tjUp'^P)> RjiTy-iip\p), Rf+tAiip' .P) 



to be solved for each submatrix R 
elements 

Ri.f^{p',p)^Ri,f^{p',p-2E,), 



JST 



composed of the 



(51) 



where 2Ep = 2-\/p2 + m? represents the collision en- 
ergy in the c.m.s. As usually, we distinguish the on- 
energy-shell {p' = p) and half-off-energy-shell {p' ^ p) 
elements. Apparently, it is pertinent to stress that even 



The first two of them satisfy the equations: 
for spin singlet 

R'jr'^{p\p)=vj't'^{p\p) 

QO 

1 '^^ T^JS=aT// „\nJS=OTt 







Ep — Eq 



Vj'r'"ip',<l)Rrj=''il,P), (55) 



^ Below in this Sect., for brevity, we shall write R(V) instead 
of RiV). 
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and for uncoupled spin triplet 



R 



and p{p) = pEp/2. 

In turn, with the help of a standard definition of the 
on-shell ^'-matrix elements 



2 J Ep — Eq 



Vft'^ip', q)Rjj='^{q.p), (56) Sf/{p) ^ exp [2*<5f^] = 1 - 2'Kip{p)Tjr ip) (60) 



while for finding the rest one needs to solve the set of 
four coupled integral equations 

Rf=,T,{p'.p) = vj^rr^,{p'^p) 

P f q^dq 



2 ./ Ep — Eg 



[vf^r}^,{p\p)Rf_tj^,{p',p) 

+VfmUp',p)RilT7-Ap\p)] 



RjiT7Mp',p) = vffr}^,{p',p) 



q^ dq 
2 ./ Ep — Eq 



mrr_,{j>',p)R^f_tl,{p'.p) 

+vjmup'^p)Rf=.f-.{p\p)\ 



Rf_tjUp',p)^vi^r}Up',P) 

P f q^ dq 



+ 



2 ./ Ep — Eq 



[vj'TF_,{p\p)Rf=,Z,{p\p) 
+vj^T}l,{p' .P)R'ji7Y+,{p\p)\ 



RiiT7+,{p',P)-vffr}Up',P) 

P f q^dq 



2 ./ Ep — Eq 



[v/^-}':,(p\p)Rf-i7Mp',p) 
-vf^r}Up\p)R'jitjUp\p)] (57) 



One should stress that after omitting the factor [1 — 
(^_iy+s+T-^ the matrix elements Vy^'^ in these equa- 
tions are given by the direct parts of our energy inde- 
pendent quasipotentials, which are derived in Appendix 
C. 

The Heitler equation (PT|) enables one to get the 
partial T— matrix elements 



■J ST It 



i7:p[p) V Ri,W{p'.P)Tilf{p.p). (58) 



I" 



whence it follows the on-shell relationship: 

Tfr{p) = RfriP) - E Rif'^ipm^riP), 



(59) 



where 



r(i?)/,f ^ ip) = T(E)/,f ^ E = 2Ep) 



the _R-niatrix elements for the uncoupled states can be 
expressed through the phase shifts 5j^: 



ST 



(61) 



Usually the isospin label is suppressed to write simply 
5 J. Of course, these quantities depend either on the 
incoming momentum p or the laboratory energy Eiab- 
We choose the second alternative. 

For the coupled states the on-shell R~ matrix ele- 
ments are conventionally parameterized in terms of the 
phase shifts 6^. and the mixing parameters ej: 



\ ^J+i J-i ^J+i J+i 
exp[2i6i] 
exp[2i6-l 



cose,/ — sine,/ 
sine,/ cose J 

cose J sine J 
— sine J cose J 



(62) 



Such a parametrization was put forward in [28]. From 
Eq.dMl) it follows that 



tan (5^ — 



^^Pip) 



R 



71 



R 



71 



and 



tan 2e j 



-"■J+1 J-1 'T- -"-J-i J+l 
f>Jl _ nJl 



J+l 7+1 ^ -"i-J-l 7-1 

R^l _i_ f/Ji 

-"•J+l J-1 ^ -"-J-1 J+l 



cose J 



(63) 



(64) 



Here in Rfi]^{p) the isospin index and the argument p 
are suppressed too. 

However, in the next section our calculations are 
shown for the so-called bar convention introduced in 
[55]. These are related to the Blatt-Biedenharn phase 
shifts by 



sin(<5^ -5i)^- 



tan 2e, 



J 



tan 2e 



J 



. n - (65) 

sin2ej ^ ' 

(cf. Ea. (|132p from [1]). It will allow us to compare our 
results directly with those by the Bonn group (in par- 
ticular, from the survey [13] )• 
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4 Results of numerical calculations and their 
discussion 

Available experience of solving integral equations sim- 
ilar to Eqs. ([55)) - ([57)) shows that it is convenient to 
employ the so-called matrix inversion method (MIM) 
[50] . [5T| (more sophisticated methods are discussed in 
the monograph |32|). In the course of our numerical 
calculations we have improved a code |33] based upon 
the MIM and successfully applied for the treatment of 
the final-state interaction in studies |34], [35] of the 
deuteron breakup by electrons and protons in the GeV 
region. 

Since we deal with the relativistic dispersion law for 
the particle energies, the well known substraction pro- 
cedure within the MIM in our case leads to equations 



Ri,f^{p',p)^v,jr{p\p) 

oo 

1 /" 
9 2^ 



2^7 P^-q" 

I- n 



{q\E, + E,)VA^,^{p',q)Ri,?,^{q,p) 



-2p'EpVi^i9,^{p\p)R'A^{p,p)] . (66) 

To facilitate comparison with some derivations and cal- 
culations from Refs. [5], [13], we introduce the notation 

= -F^ip',p) K(l', 2'; 1, 2) + Vb{2', 1'; 2, 1)] 

for the regularized UCT quasipotentials in the c.m.s. 
(see Appendix C). As in Ref [5], we put 



Fb{p\p) 



/12 



Fb[{p'-pf 



Al-{p'-pf 
Doing so, we have 

= g,^u(pO^(p) ^^;l^^), ^^ u(-pOu(-p), (67) 



and 

/ ,,/ ,,/ I „,UCT 



P f^lt^2\% |PA*lM2) = -77^^-72 — 



{p' ~p) -ml 



X jl(p') 



u{p) 



X m(-p') 



(£;p'-i?p)[7V-.9"''^ 



2m 



m(-p) 



4m^ 



{Ep, - Ep) u(p')[7o7i/ - gOfMp) 

x^i(-p')[7V-5°>(-p)^ (69) 



where (p' + p) ^ (Ep, + Ep, -{p' + p)). 



Table 1 The best-fit parameters for the two models. The third 
(fourth) column taken from Table A.l 1131 (obtained by solving 
Eas. ll66ll with a least squares fitting I36II37I to OBEP values in 
Table 2). All masses are in MeV, and n;, = 1 except for Up = 
riu, = 2. 



Meson 




Potential B 


UCT 


TT 




14.4 


14.5 






1700 


2200 






138.03 


138.03 






3 


2.8534 




An 


1500 


1200 




rrir, 


548.8 


548.8 


P 




0.9 


1.3 






1850 


1450 




fp/9p 


6.1 


5.85 




rup 


769 


769 


iO 


95/4^ 


24.5 


27 






1850 


2035.59 




■m.ui 


782.6 


782.6 


5 


9i/47r 


2.488 


1.6947 




As 


2000 


2200 




ms 


983 


983 


a,T = 


9^/47r 


18.3773 


19.4434 






2000 


1538.13 






720 


717.72 


a,T =1 




8.9437 


10.8292 




Aa 


1900 


2200 




ma 


550 


568.86 



(P»2l<f^|Pm/^2) 

.,2 



2 , , Fl{{p'~pf\ _^ 

gp,u{p )75^(p) ^^,_^^2_^2^ ^("P )75^(-p) 

(68) 



At first sight, such a regularization can be achieved 
via a simple substitution g\, — gi,Fi,{p',p) with some 
cutoff functions Fh{p' ,p) depending on the 4-momenta 
p' and p. 
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Table 2 Neutron— proton phase shifts (in degrees) for various laboratory energies (in MeV). The OBEP(OBEP*)-rows taken from 
Table 5.2 I13| (calculated by solving Eqs.JTTJ with the model parameters from the third column in Table 1). The UCT* (UCT)-rows 
calculated by solving Eas. ll66ll with the parameters from the third (fourth) column in Table 1. As in 2^, we have used the bar convention 
|29| for the phase parameters. 



State 


Potential 


25 


50 


100 


150 


200 


300 




OBEP 


50.72 


39.98 


25.19 


14.38 


5.66 


-8.18 




OBEP* 


50.71 


39.98 


25.19 


14.37 


5.66 


-8.18 




UCT* 


66.79 


53.01 


36.50 


25.27 


16.54 


3.12 




UCT 


50.03 


39.77 


25.55 


15.20 


6.92 


-6.07 




OBEP 


-7.21 


-11.15 


-16.31 


-20.21 


-23.47 


-28.70 


^Pl 


OBEP* 


-7.17 


-11.15 


-16.32 


-20.21 


-23.48 


-28.71 




UCT* 


-7.40 


-11.70 


-17.73 


-22.63 


-26.98 


-34.54 




UCT 


-7.15 


-10.95 


-15.62 


-18.90 


-21.49 


-25.41 




OBEP 


0.68 


1.58 


3.34 


4.94 


6.21 


7.49 




OBEP* 


0.68 


1.58 


3.34 


4.94 


6.21 


7.49 




UCT* 


0.68 


1.59 


3.40 


5.10 


6.52 


8.20 




UCT 


0.68 


1.56 


3.22 


4.68 


5.77 


6.68 




OBEP 


9.34 


12.24 


9.80 


4.57 


-1.02 


-11.48 


^Po 


OBEP* 


9.34 


12.24 


9.80 


4.57 


-1.02 


-11.48 




UCT* 


9.48 


12.53 


10.32 


5.27 


-0.15 


-10.27 




UCT 


9.30 


12.16 


9.81 


4.73 


-0.68 


-10.76 




OBEP 


-5.33 


-8.77 


-13.47 


-17.18 


-20.49 


-26.38 


3 Pi 


OBEP* 


-5.33 


-8.77 


-13.47 


-17.18 


-20.48 


-26.38 




UCT* 


-5.27 


-8.62 


-13.09 


-16.56 


-19.63 


-25.06 




UCT 


-5.28 


-8.58 


-12.85 


-16.06 


-18.86 


-23.79 




OBEP 


3.88 


9.29 


17.67 


22.57 


24.94 


25.36 




OBEP* 


3.89 


9.29 


17.67 


22.57 


24.94 


25.36 




UCT* 


3.86 


9.15 


17.12 


21.51 


23.47 


23.48 




UCT 


3.89 


9.25 


17.31 


21.77 


23.75 


23.61 




OBEP 


80.32 


62.16 


41.99 


28.94 


19.04 


4.07 




OBEP* 


80.31 


62.15 


41.98 


28.93 


19.03 


4.06 




UCT* 


92.30 


72.71 


51.44 


38.10 


28.20 


13.70 




UCT 


79.60 


61.53 


41.57 


28.75 


19.08 


4.60 




OBEP 


-2.99 


-6.86 


-12.98 


-17.28 


-20.28 


-23.72 


^Di 


OBEP* 


-2.99 


-6.87 


-12.99 


-17.28 


-20.29 


-23.72 




UCT* 


-2.74 


-6.43 


-12.36 


-16.54 


-19.47 


-22.78 




UCT 


-3.00 


-6.90 


-13.12 


-17.66 


-21.11 


-26.03 




OBEP 


1.76 


2.00 


2.24 


2.58 


3.03 


4.03 


ei 


OBEP* 


1.76 


2.00 


2.24 


2.58 


3.03 


4.03 




UCT* 


0.02 


-0.12 


-0.17 


0.04 


0.41 


1.40 




UCT 


1.80 


2.01 


2.19 


2.50 


2.90 


3.83 




OBEP 


2.62 


6.14 


11.73 


14.99 


16.65 


17.40 




OBEP* 


2.62 


6.14 


11.73 


14.99 


16.65 


17.39 




UCT* 


2.80 


6.61 


12.71 


16.28 


18.10 


18.91 




UCT 


2.57 


6.00 


11.32 


14.18 


15.37 


15.07 




OBEP 


0.11 


0.34 


0.77 


1.04 


1.10 


0.52 




OBEP* 


0.11 


0.34 


0.77 


1.04 


1.10 


0.52 




UCT* 


0.11 


0.34 


0.77 


1.05 


1.13 


0.64 




UCT 


0.11 


0.34 


0.75 


1.00 


1.03 


0.41 




OBEP 


-0.86 


-1.82 


-2.84 


-3.05 


-2.85 


-2.02 


£2 


OBEP* 


-0.86 


-1.82 


-2.84 


-3.05 


-2.85 


-2.02 




UCT* 


-0.87 


-1.83 


-2.82 


-2.99 


-2.75 


-1.88 




UCT 


-0.86 


-1.83 


-2.84 


-3.05 


-2.89 


-2.18 



11 





lab. energy [Mev] lab. energy [Mev] 




Fig. 3 Neutron-proton phase parameters for the uncoupled partial waves, plotted versus the nucleon kinetic energy in the lab. system. 
Dashed[solid] curves calculated with Potential B parameters (Table [TJ by solving Eqs. 1 166 II [ JTIIl l . Dotted represent the solutions of 
Eqs. II66II with UCT parameters (Table The rhombs show original OBEP results (see Table 2). 



However, a principal moment is to satisfy the re- 
quirement (86) for the Hamiltonian invariant under 
space inversion, time reversal and charge conjugation. 
A constructive consideration of the issue is given in Ap- 
pendix C. 

Replacing in equations (|67)) - (l69)) 



by 

-F^[-{p'-pr]{{p'-pr+mi}-' 

and neglecting the tensor-tensor term 

-^{Ep' - Sp)2u(p')[7o7^ - 9ou]u{p) 

x«(-p')[7V-5°'']"(-p) (70) 
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I I I I I I I I I I I I I 

50 100 150 200 250 300 

lab. energy [Mev] 

1 ' I ' I ' I ' I ' I ' r~ 




I I I I I I I I I I I I I I 

50 100 150 200 250 300 
lab. energy [Mev] 

Fig. 4 The same in Fig. 3 but for the coupled waves. 

in (|69L we obtain approximate expressions that with 
the common factor 

/Ep,Ep 

instead of 

are equivalent to Eqs.(E.21)-(E.23) from [5]. Such an 
equivalence becomes coincidence if in our formulae in- 




I I I I I I I I I I I 

50 100 150 200 250 
lab. energy [Mev] 




I I I I I 1 I I I I I 

50 100 150 200 250 
lab. energy [Mev] 



stead of the canonical two-nucleon basis |p/ii/i2) one 
uses the helicity basis as in [2]. 

In parallel, we have considered the set of equations 

^Ri,f^{p'^p) = ''v,^r{p\p) 

oo 

^p'^V,^i^,^{p\pfRt,^,^{p,p)}, (71) 
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where the superscript B refers to the partial matrix 
elements of the potential B determined in [13] with the 
just mentioned interchange of the bases. It is important 
to note that Eas. ([7T|) can be obtained from Eqs. (|55| 
ignoring some relativistic effects. In particular, it means 
that the covariant OBE propagators 



(72) 



{Ep' 



Ep? 



(P' - P)^ - ml 



are replaced by their nonrelativistic counterparts 



(P' - P)^ + ml 



(73) 



Such an approximation is a key point that gives 
rise to the potential B from T5] . 

* Sometimes associated with ignoring the so-called meson re- 
tardation (see, e.g.. Appendix E from ^2J and a discussion therein) 
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Fig. 6 The same in Fig. 5 but for the coupled waves. 



However, the transition from (1721) to (|73)) . being 
valid on the energy sheh 

Epi = Ep = -Ecms, (74) 

cannot be a priori justified when finding the R matrix 
even on the shell. Thus our calculations of the R matri- 
ces that meet the equations (1551) and (I7T]) are twofold. 
On the one hand, we will check reliability of our numer- 
ical procedure (in particular, its code). On the other 



hand, we would like to show similarities and discrepan- 
cies between our results and those by the Bonn group 
both on the energy shell and beyond it. These results 
are depicted in Figs. |3HH1 and collected in Table [51 

As seen in Figs. [SHU the most appreciable distinc- 
tions between the UCT and OBEP curves take place 
for the phase shifts with the lowest Z— values. As the 
orbital angular momentum increases the difference be- 
tween the solid and dashed curves decreases. Such fea- 
tures may be explained if one takes into account that 
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1000 
p' [Mev/c] 



1000 
p' [Mev/c] 



1000 
p' [Mev/c] 



1000 
p' [Mev/c] 



1000 
p' [Mev/c] 



1000 
p' [Mev/c] 



2000 




2000 




2000 



Fig. 7 Potentials for the uncoupled partial waves witii thie momentum po fixed as in Fig. 5. Other notations in Fig. 3. 



the approximations under consideration affect mainly 
high-momentum components of the UCT quasipoten- 
tials (their behavior at "small" distances). With the 
^-increase the influence of small distances is suppressed 
by the centrifugal barrier repulsion. 

Of course, it would be more instructive to compare 
the corresponding half-off-energy-shell i?-matrices (see 
definition (jSip ). Their p'-dependencies shown in Fig.5- 
6 are necessary to know when calculating the ■(/j^^^ scat- 
tering states for a two-nucleon system. Such states may 



be expressed through the partial-wave functions ip^, a = 
J, S, T that have the asymptotic of standing waves (see, 
for example, [34]). Within the MIM every (p can be rep- 
resented as 



N+l 



where the coefficients B^,{j) are the solutions of the 
set of linear algebraic equations approximately equiva- 
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Fig. 8 The same as in Fig. 7 for the coupled waves. 



lent to the integral equation for the corresponding R- 
matrix; A'^ is the dimension of the set, pj are the grid 
points associated with the Gaussian nodes on the in- 
terval [—1, 1], pn+i — Po (details can be found in [Ml 
I35j). IVIeanwhile, our computations presented here have 
been done with iV = 32 (we do not talk about tests with 
other A^- values to get the results stable with respect to 
N). 

Completing the comparison of the UCT and Bonn 
results we display in Figs. 7-8 some interactions in rel- 



evant partial states. Appreciable distinctions between 
dotted and solid curves in these figures mean that the 
UCT and Bonn parameters from Table 1, ensuring a 
fair treatment of such on-energy-shell quantities as the 
phase shifts, may be inadequate in constructing model 
nucleon-nucleon potentials. It seems to be especially 
prominent in case of the '^Si Di potentials respon- 
sible for the formation of the tensor part of nuclear 
forces. 
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In addition, one should emphasize that hitherto we 
have explored the OBEP and UCT i?-matrices in the 
c.m.s., where the both approaches yield the most close 
results. It is not the case in those situations when the 
c.m.s. cannot be referred to everywhere (e.g., in the 
reactions NN — )■ 7 A'' A'' and -yd pn). In this respect 
our studies of the differences between UCT and OBE 
approaches are under way. 



5 Summary, Conclusions and Prospects 

The present work has been made to develop a consis- 
tent field-theoretical approach in the theory of nucleon- 
nucleon scattering. It has been shown that the method 
of UCT's, based upon the notion of clothed particles, 
is proved to be appropriate in achieving this purpose. 

Starting from a primary Lagrangian for interacting 
meson and nucleon fields, we come to the corresponding 
Hamiltonian whose interaction part Kj consists of new 
relativistic interactions responsible for physical (not vir- 
tual) processes in the system of the bosons (tt— , ry— , p— , 
w— , 5— and cr— mesons) and the nucleon. Proceeding 
with the CPR we have confined ourselves to construct- 
ing the four-legs interaction operators K\ ' in the two- 
nucleon sector of the Hilbert space of hadronic states. 
The corresponding quasipotentials (these essentially non- 
local objects) for binary processes ISSISS — >■ NN ^ NN — >■ 
NN, etc. are Hermitian and energy-independent. It 
makes them attractive for various applications in nu- 
clear physics. They embody the off-shell effects in a 
natural way without addressing to any off-shell extrap- 
olations of the 5— matrix for the NN scattering. 

Using the unitary equivalence of the CPR to the 

BPR, we have seen how in the approximation Kj = 

Kj ' the extremely complicated scattering problem in 

QFT can be reduced to the three -dimensional LS- 

type equation for the T— matrix in momentum space. 

The; ciquation kernel is given by the clothed two-nucleon 

interaction of the class [2.2]. Such a conversion becomes 

(2) 

possible owing to the property oi Kj to leave the two- 
nucleon sector and its separate subsectors to be invari- 
ant. 

Special attention has been paid to the elimination 
of auxiliary field components. We encounter such a ne- 
cessity for interacting vector and fermion fields when 
in accordance with the canonical formalism the inter- 
action Hamiltonian density embodies not only a scalar 
contribution but nonscalar terms too. It has proved (at 
least, for the primary pN and loN couplings) that the 
UCT method allows us to remove such noncovariant 
terms directly in the Hamiltonian. To what extent this 



resiilt will take place in higher orders in coupling con- 
stants it will be a subject of further explorations. 

Being concerned with constructing the two-nucleon 
states from H and their angular-momentum decompo- 
sition we have not used the so-called separable ansatz, 
where every such state is a direct product of the cor- 
responding one- nucleon (particle) states. The clothed 
two-nucleon partial waves have been built up as com- 
mon eigenstates of the field total angular-momentum 
generator and its polarization (fermionic) part expressed 
through the clothed creation/destruction operators and 
their derivatives in momentum space. 

We have not tried to attain a global treatment of 
modern precision data. But a fair agreement with the 
earlier analysis by the Bonn group makes sure that our 
approach may be useful for a more advanced analysis. 
In the context, to have a more convincing argumenta- 
tion one needs to do at least the following. First, show 
the low-energy scattering parameters and the deuteron 
wave function calculated within the UCT method. Sec- 
ond, consider triple commutators [R, [R, [R, Vb]]] to ex- 
tract the two-boson-two-nucleon interaction operators 
of the same class [2.2] in the fourth order in the coupling 
constants. Third, extend our approach for describing 
the NN scattering above the pion production thresh- 
old. All the things are in progress. 

As a whole, the persistent clouds of virtual particles 
are no longer explicitly contained in the CPR, and their 
influence is included in the properties of the clothed 
particles (these quasiparticles of the UCT method). In 
addition, we would like to stress that the problem of the 
mass and vertex renormalizations is intimately inter- 
woven with constructing the interactions between the 
clothed nucleons. The renormalized quantities are cal- 
culated step by step in the course of the clothing proce- 
dure unlike some approaches, where they are introduced 
by "hands". 
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A Model Lagrangians and Hamiltonians 
within Canonical Formalism. Transition to 
Clothed Particle Representation 

We will focused upon three boson fields (pseudoscalar, scalar and 
vector) that are coupled with the nucleon by means of the often 
employed interaction Lagrangians, 

£s = -g°^^^s, (75) 
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4m 



(76) 



(77) 



i(7M7^ - 7.^7^) and = 9"^? - S'^- As 



where ct^ 

Refs. f 1101 . nil ), throughout this paper we use the definitions 
and notations of [25], so, e.g., 7^7^ + 7^7^^ = 2g^i, {goo = 1, 

gil = 922 = 933 = -1), 7^ = 707m 70 , 75 = i7°7^7^7^- In addi- 
tion, following Ref. I12| we shall distinguish via upper(lower) case 
letters between the Heisenberg and Dirac picture field operators 
(ip" and <P vs ip and ip, respectively, for fcrmions and bosons). 

Of course, we could incorporate the so-called pscudovcctor 
(pv) coupling 



Ulps 



(78) 



Since in this paper all used model interactions are suggested to be 
modified by introducing some cutoff factors, it has no matter that 
couplings II77I I and II78II with derivatives are nonrenormalizable 
(cf. an instructive discussion of this subject in Subsect 3.4 of 
the survey |13l V In the context, starting from couplings II75II - 
II77I I with "bare" constants Qp^, g^, g^ and we have tried 
to reproduce some results obtained in Refs. [2], [5], where such 
constants from the beginning are replaced by effective parameters 
9ps, <?s, 9v and /v. It explains (at least, ad hoc) our restriction to 
these Lagrangian densities. Recall also that for isospin 1 bosons 
one needs to write instead of ^j,, where t is the Pauli vector 
in isospin space. 

In constructing the Hamiltonians with Lagrangian densities 
||75| |- ||77|| as a departure point, we have first used the equations 
of motion for the H fields and the so-called Legendre transforma- 
tion (from L to H) to express the total Hamiltonian H in terms of 
the independent canonical variables and their conjugates. Then, 
passing to the D picture (interaction representation) the Hamil- 
tonian has been split into a physically satisfactory free-field part 
Hq and an interaction V. Int. al., since the component has 
no canonical conjugate, we have resorted to a trick prompted by 
Eq.(7.5.22) from |12| to introduce a proper component c/jJ in the 
Dirac picture. As a result, we arrive to the interaction Hamilto- 
nian densities: 



of the interaction density 'K(x) in the D picture 

■K{x) = 'Ks{x) + -Kps{x) + 'K^{x), (85) 

taken at t = 0, i.e., in the Schrodinger (S) picture. In other words, 

:k(x) = :k(o,x). 

Expressions Il79l l-| [84ll exemplify that for a Lorentz-invariant 
Lagrangian it is not necessarily to have "... the interaction Hamil- 
tonian as the integral over space of a scalar interaction density; 
we also need to add non-scalar terms to the interaction density 
..." (quoted from p. 292 of Ref. I12I V It is the case with derivative 
couplings and/or spin > 1. 

In this connection, let us recall the property of the density 
'K{x) to be scalar, viz.. 



U f{j\, a)'K{x)Up'^ {A, a) = H(Ax + a), 



(86) 



where the operators Up{A,a) realize a unitary irreducible rep- 
resentation of the Poincare group in the Hilbert space of states 
for free (non-interacting) fields. Other comments will be given in 
Appendix B. 

Moreover, it is well known that this property is a key point of 
covariant perturbation theory of the 5— matrix (see, e.g.. Chapter 
V of 1121 and refs. therein). In the respect, the division JSTJ is 
not accidental. 

Further, the BPR form ([1} and its comparison with Eq. II84I I 
give rise to the interactions itSj-© between the bare bosons and 
fermions with physical masses. The next step is to apply the 
clothing procedure exposed in I10| (see also 1111 ). where the first 
clothing transformation VF'^' = exp[i?(-'^^] (ij'-*^' = — i?'-'') elim- 
inates all interactions linear in the coupling constants. Its gener- 
ator i?'^' obeys the equation 



with 

=Vs + Vps + vy\ 

where 



(87) 



■K,{x) = g'ii,(x)^(x)^s{x), 



(79) 



V'(x)7mV'(xX(x) 

/v 



-Hf^V(x)^M-i/'(x)<prW . (88) 
Am ' 



"yipsix) = igp^ip{x)'ysip{x)ipps{x), 
'K^(x) = -K^a^{x) + -K 

ncov 

where 



(80) 
(81) 



J{cov(x) = <;Xa;)7^i^(2:)vp(;(x) + ^,/.(a;)a^,V(2^)¥'r {^) (§2) 



f° 

J V 

4m 



and 



Mncov(a;) = -^-^tp{x)-^o^{x)ip{x)'^otp{x) 
/02 _ 

+ --^^{x)(Toi^{x)^{x)croi^(x). (83) 
4m^ 

It is implied that the total Hamiltonian of interest H = Ho + V 

consists of the sum Ho = Hoj + ^0,6: where Hqj the free- 

b 

fermion Hamiltonian, Ho b separate free-boson contribution, and 
the space integral 



X) 



(84) 



Following [TD] Eo. llHTt is satisfied with 



_r(i) = -i lim 
e->a+ 



00 

j Vlj\t)e-^*dt 



(89) 



if m;, < 2m. 

The corresponding interaction operator in the CPR (see re- 
lation II13I I of the text) can be written as 



i^f)^ 1 ^^(l)^^(l)j +y(2)_ 



(90) 



where we have kept only the contributions of the second order in 
coupling constants, so 

1/(2) ^ j |^^V;(x)7oi/'(x)?/;(x)7oi/'(x) 

-f^V;(x)ao«^/'(x)i/;(x)f7o,i/'(x)| (91) 

that coincides with the space integral of the non-scalar density 
'Jincov{x) at t = 0. It is implied that all operators in the r.h.s. 
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of Eg. 1190 II depend on the clothed-particle creation(destruction) 
operators etc = W^^^^ aW^^\ for example, involved in the Fourier 
expansions 



V'(x) = (2^)-3/2 / dp f^J2 HPf^)bc{pt^) 

J V ^ 

+t;{-p/i)4(-PM) exp(ipx) 



+e''(-k, s)aj(-k, s)j exp(jkx), 



where the e''(k,s) for s = +1,0,-1 are three independent vec- 
tors, being transverse fc^e^'(k, s) = 0, and normalized so that 



^ e^(k, s)e* (k, s) = -g^^ + kf_ikv/r 



In this paper we do not intend to derive all interactions between 
the clothed mesons and nucleons, allowed by formula KOII . Our 
aim is more humble, viz., to find in the r.h.s. of Eg. I I 901 1 terms 
of the type II18II . responsible for the N-N interaction. Along the 
guideline the commutators ^ Vs] and \ ,Vps\ generate 

the scalar- and pseudoscalar-meson contributions Ki,{NN — > 
A'^A'^) with coefficients II20II - II21II . In case of the vector mesons 
we encounter an interplay between the commutator i[iJv'^\ K/^'] 
and the integral II91II . 

To show it explicitly, let us write 



(1) 



-L_5:/dkdp>- 



-.e''{k,s) 



(27r)3/2 ^2uj^E^,Ep 
(5 (p' - p - k) m(p'm') |gv7p - ^icTvpk"^ u{pfj.) 

X bl{p'fi')bc{pfi)ac{K s) + H.c, (92) 



with 



K^(NN -s> TV A) 



^(^E j dp\dp'^dp^dp2 
& (Pi + P2 - Pi ~ pa) 



(pi - Pi)^ - ml 

/v 



"(PiM'i) "i (Sv + - ^(P'l \ "(Pl/^l) 



X uipl^pi'^) |(gv + /v)7" - ^{P2 +P2)"'| «{p2/^2) 



P2 



«{PiVi) |(Sv + /v)7i' - ^{P'l +Pi)i'| "(PiMi 



x«(p^/.^)A{7V 
2m 



tt(p2M2) 



X ^'J(Pl/i'l)f'J{P2/^2)''c{Pl/^l)^c{p2/i2) 



and 

ATcontCATA ^ AAf) 



/-^p'l^P2^Pl-^P2 



^(^E j dp'^dp'.,dpidp2 

^ (Pl +P2 - Pi - P2) 



^"{PlA'l)70«(Pl/il)«(P2A'2)70«(P2/i2) 



Am? 



"(PlVi)7o7"(PlMi)«(P2M2)707«(P2M2) 
X ^'J(PlVl)feJ{P2/^2)^c(pi/il)bc{p2A»2) 



The latter may be associated with a contact interaction since it 
does not contain any propagators (cf. the approach by the Osaka 
group 1391 ). It is easily seen that this operator cancels completely 
the non-scalar operator V^^\ In other words the first UCT en- 
ables us to remove the non-invariant terms directly in the Hamil- 
tonian. It gives an opportunity to work with the Lorentz scalar 
interaction only (at least, in the second order in the coupling con- 
stants). In our opinion, such a cancellation, first discussed here, 
is a pleasant feature of the OPR. 

The remaining vector— meson contribution Kv{NN —> NN) 
is determined by coefficients II22II and gives us one more relativis- 
tic interaction in the boson-fermion system under consideration. 



(2^)3/2/^7 ^V2(^ki-p'A'p 

5(p^-p-k) r /v. \ 

X ^ "(P fJ- ) i g^lp - —i(Tupk \ u(p^i) 

-Bp' - -Bp - ojk I, 2m J 

X bl{p' ^l')bc{p^l)ac{]^, s) - H.c, (93) 

retaining only those parts of y(i) and R(^\ which are necessary 
for deriving K^^NN — NN). After a simple algebra we find 

- [-R(i',y(^)l (AA-!> AAf) 
2 L -I V 

= K^{NN NN) + K^ontiNN NN) 



B Partial Wave Expansion of Two Particles 
States in the CPR 

We will show how one can proceed without the separable ansatz 
|PlP2MlM2> = IpiMi) |P2/-t2> 

often exploited in relativistic quantum mechanics (RQM) (see, 
e.g., |23l and [241 ) in getting expansions similar to Eq.(47). Unlike 
this, our consideration with particle creation/destruction as a 
milestone, where the clothed two-nucleon state is given by 

|Pl/'i;P2/*2> = 6l(piA'l)4(P2M2) \^) (94) 

^ Sometimes it is convenient to handle the operators bc{pfi) = 
^/pofeJ(p/i) and their adjoint ones bc{pfJ,) that meet covariant re- 
lations hl{p'fi'),ba{pfJ.)\ =po'5(p' - p)Vm 
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By definition, it belongs to the two-nucleon sector of IK being the 
Kp eigenstate with energy i? = +p§. Moreover, it is assumed 
that vector being the single clothed no— particle state, has 

the property 



UF(A,a) \n) = 



\n) 



(95) 



V yl G L+ and arbitrary spacetime shifts a = (a'^ , a) 

to be invariant with regard to the Poincare group space 
inversion and other symmetries. Here is the homogeneous 
(proper) orthochronous Lorentz group. In turn, every Up{A, a) is 
expressed through the clothed free-particle generators of space- 
time translations 

Pl^{a,) = {Kp{a,),FF{a,)}, 

space rotations 

Jpia^) = {M|.3(a,),M|.i(ac),MFK)} 

and the Lorentz boosts 

Bpia^) = {A/0i(ac):^^F^K),Af03(Q,^)} 



Up (j^, o.) = exp 



(96) 



where the antisymmetric tensor uj^u = ^i^vi^ with {fj,, = 0, 1, 2, 3) 
has six independent components. 

Now, let us remind of the following transformation property: 

Up {A, a) fet (pf,) Up' {A, a) = e'^^-^bl {Ap,,') D^^/^^ (W (A,p)) 

(97) 

with the D'^/^) function whose argument is the Wigner rotation 
W(A,p). The latter has the property W{R,p) = R for any three- 
dimensional rotation R. In other words, from Eq. |I97I I it follows 
that under such a rotation, when Up{A, a) = Up{R, 0) = Up{R), 
one has 



Up (R) bl (pM) [7-1 (R) = bl {Rpp,') dIV^^ (R) 



(98) 



In addition, we need to have an analytic expression for the op- 
erator 3 f firm to be expressed in terms of b], and be- To do it we 
recur to the well known result: 



' ferm 



with S : 



I 



diiip* (x) 



d_ 



1. 



^(x) 



(99) 



(T 
<T 



where <t the Pauli vector. 



After modest effort one can see that 



' ferm 



-'ferm 



^ ferm i 



(100) 



where lif^rm (S/erm) the orbital (spin) momentum of the fermion 
field, that are given by 



-'ferm 



j dppx — ^- — 6c(pm) -6I(pm) — — 



^'''(p'^)d.(pM)-4{PM)'^''^^''^^ 



dp 



dp 



(101) 



® The correspondence (A, a) — > Up{A,a) between elements 
{A, a) £ n and unitary transformations Up{A,a) realizes an ir- 
reducible representation of i7 on J{ (see, e.g.. Chapter 2 in |12p 



and 

^ ferm ' 



dp ■ 



u\p,^i')'Su(p,^i) 



ip X u^p,^^') 



, du{p,fi) 9ut(p,p') 



dp 



dp 



"(P./»)j I ^l(P/i')^c(pM) 



vHp,^,')^vip,^,)~^px (vHp,t^')^^^^ 



'J(p>m) I \ dl{pfi')de{pfJ.) 



dp 
or 

S/erm = \ ^ j '^'P x],/<TXm { ^1 (PM')*c (p/i) - dl{pp.')dc{p^)^ , 

(102) 



t/ ,s<9"(P:M) 3Mt(p,;i') 

u'(p,^l ) — m(p,m) 



dp 



dp 



" (P: M ) ^ H '^(P. W 

op op 



xL^Xm X P. 



if one uses the Dirac spinors defined in "25; . In these formulae we 
sec the Pauli spinors x^l■ It is important to keep in mind that 

and 



iferm 1^ \ j d^X^t (x) S,/.(x) 



In fact, the operator Sf^^rm stems from a destructive interplay 
between orbital and spin parts of decomposition II99I I. In our 
opinion, such an interpretation differs from the definition below 
Eq.(13.48) in 

Hitherto, we have preserved the Oc-dependence of the rele- 
vant operators although, as shown in Sect. 3 of 1101 . in the instant 
form of relativistic dynamics the operator F p{ae){3 pictc)) co- 
incides with the total linear momentum P = 'Pp{a) (the total 
angular momentum J = 3p{a)). Such an observation allows us 
to omit the label c if it does not lead to confusion. 

After these preliminaries we prefer to proceed sufficiently 
straightforwardly repeating the well known steps (cf. |23| and 
refs. therein). First, when handling the c.m.s. two— nucleon state 




|PM1M2> = Ipmi; -PP2> 

let us consider the vector 



\pSMs) = ( \i^l\t^2 



|PMlM2> = ( ^Ml^M2 



SMs |PM1M2> : 



SMs \pSMs) 



(103) 



(104) 



(105) 



^ Of course, what follows can be extended to an arbitrary 
frame. 

* For a moment, the isospin quantum numbers are suppressed. 
We will come back to the point later. 
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Second, one introduces 

\pJ{lS)Mj) = J dpYi^^ (p) \pSMs) {ImiSMs \JMj ) (106) 
or, reversely, 

\pSMs) = [ImiSMs \JMj ) Y*^^ (p) \vJ{lS)Mj) (107) 

with the unit vector p = p/p. 

Third, substituting II107II into the r.h.s. of Ea. lll05ll we arrive 
to the desired expansion, 



C The Regularized Quasipotentials and their 
Angular Momentum Decomposition 

Trying to overcome ultraviolet divergences inherent in solving 
equation II50I I. we will regularize their driving terms by introduc- 
ing some cutoff factors. It can be achieved if instead of Eg. 1191 1 
one assumes 



V-req / / / / / \ 
6 (Pl/il.P2M2;PlMl:P2/i2) 



^ (Pi + P2 - Pi - pa) 

VP01P02POIPO2 



|PM1M2> 



1 1 
-M1-M2 



SMs {ImiSMs \JM.,)Yi' (p) \pJ{lS)M.j) 



(108) 



From the physical viewpoint it is important to know that 

J 2 \p.J{lS)Mj) = J{.J + 1) \pJ(lS)Mj) (109) 

\pJ{lS)Mj) = Mj \pJ{lS)Mj) 
and 



f erm 



f erm 



\pSMs)=S{S + l)\pSMs) 
\pSMs) = Ms \pSMs) 



(110) 



A simple way of deriving these relations is to use the transforma- 
tion 



UpmpSMs) = \RpSM's)D(^}^^^^^iR} 

V -R S the rotation group 
and its consequence 

Uf{K)\pJ(IS)Mj) = \pJ(lS)M'j)D['^]^j^_^^iR) 



(111) 



(112) 



for infinitisemal rotations. With the help of lllOI I it is easily seen 
that 



3 \pJ{lS)Mj) = S{S + 1) \pJ{lS)Mj) 



(113) 



Thus we have built up in the CPR the common eigenvectors of 
operators and S^^^^. Probably, one should note that S is not 
any eigenvalue of the so— called invariant (or internal) spin oper- 
ator I = (0,1) introduced by M. Shirokov [26|. While the latter 
involves internal orbital motion and polarization contributions, 
the quantum number S, whose values are regulated by Clebsch- 
Gordan coefficient in definition II104II , characterizes rather the to- 
tal spin of the two-nucleon system. Coming to the end, we allow 
ourselves to write the parity operator of the fermion (nucleon) 
field in the CPR (cf. Eq.(15.93) from [25]): 

■^/erm ^■^P['^'^ ferm]^ 

■^ferm = ~\ j '^^P [^1 (P/^)''c (P/i) " (p/i)fec (- P/i) 

+dUpp)dciptJ.) + 4(pM)dc(-pM)] (114) 

with the relations 

y/e™4(P/^)5'7e™ = -4{-PM) 

that extend the rules II97II for the clothed particle operators to 
such an improper transformation as the space inversion. 



omitting for the moment isospin indices, so 

Kt{NN -s> NN) -5> Kl'^'^iNN -s> NN) 

[,,,,, ^dpl dp,^ dpi dp2 

= s [Pi + P2- Pi - P2) -r--r- 

J Poi P02 POl P02 

(PlA'l.P2/^2;Pl/^l:P2At2) 



fej (P>'l ) (P2 M2 )fec (Pl Ml )fec (P2 M2 ) (116) 



(117) 



where the old ones vi, are determined by Ea. ll20ll - ll22ll and empiri- 
cal cutoff functions should not violate the known symmetries 
of interactions. In particular, if one writes 



(NN -> NN) = J Kl^^{-x)dn 



(118) 



with 



{27tY J p,j2 POl P02 

ffc (PlMl>P2M2;PlMl:P2M2) 

f'c(piMi)fec(P2M2)f'J(Pl/^l)^'J(P2M2) (119) 
then the RI implies the property of the operator 
Kl^'^x) = exp {iKpt) Kl^^yi) exp (-iKpt) (120) 
to be a scalar, viz., 

Uf(A,X)KI''<>{x)U-\A,\) = K^^i^iAx + X) (121) 
But accordingly I I97II it imposes the following restrictions 



Dlil,{W(A,p'i))Dlil, (WiAy^)) 



X Di,}> iW(A,pi)) D^ll; {WiA,p2)} 



X'"b (PlMl:P2M2;PlMl>P2M2) 
r 



vl''^iAp'i7j[,Ap'2V2;Apir,i,Ap2V2) (122) 



needs to verify that the old (non-regularized) ones satisfy relation 
112211 themselves. It can be done with the help of the property 

S{A)u{pi.) = DlV^^\w{A,p))u(Ap^,'), 

^ We do not consider the functions depending on nucleon po- 
larizations 
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where S{A) = exp [— ^oJ^i/cr^"] the matrix of the nonunitary 
representation A S{A) in the space of spinor indices. Indeed, 
recalhng more the relations 

S{Ar'YS{A) = 'y''Af^, 

one can easily seen that the quantities vi, (p j ^'j^,p2M2iPlMl'P2M2) = 
v,,{l',2'; 1, 2) by Eqs.(20)-(22) obey Eg. lIT^ . Also, let us remind 
(see, e.g., |27p that for a Lorentz boost A = L{v) with the veloc- 
ity V, the Wigner transformation W{A,p) is the rotation about 
the V X p -direction by an angle -0, which can be represented as 



Po 
m 



PO 
m 



: 7|V X p| , 



where the zeroth component of the nucleon momentum p* = 
(pg,p*) = L{v)p in the moving frame and 7 the corresponding 
Lorentz factor. As noted, W{R,p) = R for a pure rotation R. 

Now, keeping in mind the relation I I117I I we need to deal with 
a Lorentz-invariant cutoff, 



Fb{p'i,p'2;Pi,P2) = Fi,{Ap[,Ap'2;Api,Ap2) 
in our model regularization. 



,.reg, , , , , ^ (Pl' + P2 " Pi " P2 



(123) 



VP01P02POIPO2 
-P'i, (P'l . P2 ; Pi . P2 (P'l /i'l , P2 At2 ; Pi Ml : P2 /i2 ) (124) 

Further, assuming that Fj, depends on the two invariants (pi — 
p'j2 and (p2 - p^)2, 

Ftip[,p'2;pi,P2) = h ((Pi - Pi)^ (P2 - P2)') 
we have on the momentum shell, 
(pi-p'l)2 = (Si-i?()2-(pi-p^)2 
^P2-p'2? = (E2-E!,f-(p^-pif 

SO in the c.m.s. one has to deal with the function /(,[(p — p')^] 
(^6 [(p' — p)^] in the main text) of one invariant (p — p')^ = 

(Bp-EpO'-(p-p')'- 

Of course, a similar regularization could be implemented if 
the vertices 



the corresponding prediction for the FF Gt^nn (e.g., in an an- 
alytic form put forward in 1401 ) could be employed as the bare 
factor g^j^{p' ,p)- The latter has the property g^^{Ap' , Ap) = 
g^j_{^' ,p), i.e., it should be dependent upon the Lorentz scalar 
p'p or (p' — p)^ to fulfill condition (86). Other vertex cutoffs g^_, 
g-+ s-nd 9-_ ^I's also desirable to be introduced on the same 
physical footing. 

By the way, along with properties g'^_j_{Ap' , Ap) = g^-j- (p' ,p) 

condition (86) means that the crossed cutoffs g^j {p',p) and 

g^^(p',p) should be some functions of the four-product p'_p = 
p'p—. In turn, other symmetries, mentioned below Eq.(69), yield 



the following links 3°, ^(p',p) = g^,^(p,p'), 9++(p',P) - 



.(p',p), 



g^j {p',p) = g^_^{p,p') for the real functions. All the things are 

extremely important for constructing a relativistic nonlocal QFT, 
where the boosts operators are determined as elements of the Lie 
algebra of the Poincare group. It will be presented somewhere 
else. 

In addition, let us address the matrix elements (p' | Ji,(0) | p) 
that are closely connected with the FFs in question (see, e.g., 1141 ) 
since we handle the corresponding baryon current density Ji,(x) 
at X = sandwiched between physical (clothed) one-nucleon 
states. Such matrix elements might be evaluated in terms of the 
cutoffs and other physical inputs using some idea from I41| (cf. 
the clothed particle representation of a current operator therein) . 
From the constructive point of view, it means the use of defini- 
tions I p) = 6c(pm) I ^} a-nd I p'} = 6|(p'/i') I 17), adopted at the 
beginning of Appendix B, and the similarity transformation 

= e«("<=)j-(0)e-«("'^) 
= ig(0) + [R{a,),jim + i [i?(ae), [i?(ac),ig(0)]] 

-f I [/?(«.), [iJ(a,),[iJ(ac), jX(0)]]] + ... , 

where i^(0) is the same current density but expressed through the 
clothed operators. The nonperturbative expansion in the commu- 
tators gives an opportunity for a systematic evaluation of correc- 
tions to matrix elements 

{n I b,ji{<S)bl I Q) = {Qo 1 6i6(o)fet I no). 



ii-{p'ij-')rb{k)u{pii) 

(see, e.g., Eq.(92)) would be at the beginning modified by intro- 
ducing cutoff factors g\^{p' ,p) with the bosons and fermions on 
their mass shells = m^, p'^ = p^ = and the momentum 
conservation k = p' — p. Neglecting a possible dependence of such 
factors on particle polarizations, we preserve the subscript H — h 
keeping in mind a more symmetrical consideration, where some 
factors g^,^{p',p) could be introduced with the proper energy- 
sign labels e' and e for separate three-legs contributions to a 
given boson-fermion Yukawa-type interaction. Then the corre- 
sponding regularized quasipotcntials would be contained factors 



-F'b(p'l.P2;Pl>P2) = 5++(pi,Pl)3++(p2,P2) 



(125) 



in the quadrilinear terms. 

As noted below Eq.(66) our choice of the phenomenologi- 
cal cutoff -Fi,[(p' — p)'^] has been prompted by those investiga- 
tions [2] of the Bonn group. Of course, recent developments in 
studying the structure of meson— baryon interaction vertices are 
of great interest. Among certain achievements in this area we 
find a microscopic derivation of the strong wNN and ttNA FFs 
within the quark model developed by the Graz group (see, e.g., 
|40|). Being free of any phenomenological input (fit parameters) 



Some simplifications originate from the well— known fact that sim- 
ilar expectations of the commutators that involve odd number of 
meson operators are equal to zero. 

In general, one can elaborate a recursive procedure of cal- 
culations, like that by Kharkov-Padova group, for manipulations 
with the multiple commutators [V]" (n = 2, 3, • • • ) (see |llp. Do- 
ing so, one can find corrections to formula (125) obtained from 
the commutator [V]-^ = [R, V]. This work is in progress. 

After this prelude we note that the partial-wave matrix ele- 
ments of interest are defined by 

b b 

{I'm'iSM'g \JMj) {ImiSMs \JMj ) {p'SM'g\ Vt \pSMs) 

(126) 



with 



{p'SM's\Vt\pSMs) 



Vi,{p'SM's,pSMs) - V,',,{p'SM's,pSMs), (127) 
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2 2 



-2 [(Bp, + £;p)G2 + G3 + (Bp, - Ep)Ga\ } , (137) 



(128) 



where we have employed formula (45) and property (48). 



In turn, the matrices ?)*''(p'/ij/X2, p^i/J2) can be represented as 



= ri(pVi;pMi)D6(p',p)n(-pV2;-PM2), 



Db(p ,P) = —, 2 ■ 

(P - P) - ml 

Recall that 



so 



(129) 
(130) 



V^^,(p'SM's,Y>SMs) = {-lT + 'V!„{p' SM's,-pSMs). (131) 

Further, taking into account the completeness of the matrices 1, cr 
in 2 X 2 space, the non-regularized vertices from Eas. ll20l l- ll22l l 
can be written as 

Bj,(p'/i'; p/i) = Ai,(p '; p)<5j,, J, + a hnoar functional of cr 

Now, to get the matrix ( p ' SM'g | Vi, IpSAfg ) we could do all sum- 
mations in formula II128II over ^ projections directly. However, we 
prefer the following way putting formally 

(/iklMi) = (MI<t(1)Imi> and (/i2|(T|^2> = (At2l<^(2)|M2> 
to obtain 



(132) 

with the SMs eigenvalue equations, 
S^\SMs) = S{S+l)\SMs) 
SalSMs) = Ms\SMs), 

where S = ^ [^(1) + ''■(2)]. The operators Gj can be expressed 
through the operator S with the help of the relations 

<t(1) • n <t(2) ■ n = 2 ■ n)^ - 

o-(l) X n <t(2) X n = ■ <t(2) - <t(1) • n a-(2) ■ n 
for any vector n. 

As a result, we find with the models II20I I- II22I I 



Gs 



rCDsip',p) {U^+Vl - 2(V3S) [iUi + (VaS)]} , 



G 



- ^^^^"'^^''^^ " 2(ViS)2} , 



f~lVV I /^Vt I /^tt 



9l 



2(27r)' 



rCDv(p',p)Gi 



(134) 
(135) 

(136) 



G: 



1 



CD^{p',p) {4m2Gi 



2(27r)3 4m2 
- 4m [{Ep, + Ep)G2 + Ga + (Bp, - Bp)G4] 

+G5 + 2(Bp, -Bp)G6}, (138) 

with 

Gi = Ul + vi + - 2(8^ - l)Vi 

+ 2(V3S) [iU2 - (VaS)] + 2(ViS)2 + 2j(V2 x Vi)S, 

G2 = U1U2 - vi + 2(V3S) [i(Ui - U2) + 2(V3S)] , 

G3 = C/1BV2 + j(BV2)(V3S) + i(7i(B X Vi)S 

- (B X Vi) [2S(V3S) - V3 - i(V3 X S)] , 

G4 = (5 - 2S2)ViV2 + (ViS)(V2S) + (V2S)(ViS), 
G5 = b2 {Ul+Vi + 2(V3S) [if/i - (V3S)]} , 

Ge = C/iBVi + j(BVi)(V3S) + iC/i(B x V2)S 

- (B X V2) [2S(V3S) - V3 - i(V3 X S)] , 

where we use the notations 

(7i(p',p) = (Bp, -(-m)(Bp +m) -p'p = Ui, 

'^2(p',p) = {Ep, + m){Ep + 'm) + p'p = U2, 
Vi(p',p) = p(Bp, + m) - p'(Bp, +m) = Vi, 
V2(p',p) =p(Bp, + m) + p'(Bp + m) = V2, 
V3(p',p) = p'xp = V3, B(p',p) = p'+p = B 
and 

C = [4Bp,Bp, (Bp, + m)(Bp + m)] "\ 

These expressions were used by us to evaluate the matrix ele- 
ments of interest. 



X ( yjiij (P ')|G6(P ', P; S) + (-l)^G5(p ', -p; S)| Y'ji,^ (p)) 



2J- 



X (F5'S,(p')|G6(p',p;S)|F5i,^(p)), (139) 



where 

^IS 



(^33) I Y'jI,^ (n)> = y,^, (h) {l'm[SM's \JMj ) \SMs) 



the so-called spin-angular states. 

A simple extension to the states with the isospin T yields 
the factor h{T) [l - (-1)^+'+^] instead of [l -|- (-1)'5+'] in the 
r.h.s. of Ea. lll39ll . Its appearance results in the well— known selec- 
tion rule for the nucleon-nucleon scattering. Here 



1 for neutral bosons, 

2T{T + 1) — 3 for charged bosons 



(140) 
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The operators G[,(p',p;S) depend on the scalars (p'p), (Sp')^, 
(Sp)^, (Sp')(Sp), (Sp)(Sp'), so all we need is to calculate the 
following integrals: 

= J dp'|dp(y5'4(p')|i?b(p',p)/{p',p;S)|yji^^(p)>, 

where I{p' , p; S) is a polynomial of these scalars. 

After a lengthy calculation we arrive to our working formulae. 
In particular, in case of the tensor-tensor interaction in the p- 
exchange channel we have for uncoupled waves 

{[pV + (P^ +p'^)(E'E - 2m2) + 6m^{E'E - m^)] PQj{p',p) 
-pp'(p2 + p"" + 4m^)PQ^j\p' ,p) - p''p'^PQ'-f\p' ,p) 



-(E' - Ef [{E'E - 5m^rQj{p',p)-pp'PQ'-j'>(j}',p)] } . 



where P„(cos6) is the Legendrc polynomial. Using the Neumann 
integral representation for the Legendre function of second kind 
Qnix) one can write for any n;, 



Qi(p',p) = 2n J d(cose)— ^ 



(cos0) 



: — 47r 



[4- 



(2p'p) 



(p' - pY - ml 

1 \{y-~ 



2ni, 



2nfc-l 



My) 



y = 

and 



2p'p 



2p'p 



(1^1) QUoip',P)^0. 
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